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Abstract 

Although the problem of tire modeling and analysis has been a subject of 
continuing concern for the tire industry, to date no simple and general tire model 
exists for predicting the response of the tire under various loading conditions. 

Much of the recent progress in finite element technology has not been exploited for 
tire modeling and analysis. The present paper focuses on this issue. Specifically, 
the paper reviews some of the recent advances in finite element technology which 
have high potential for application to tire modeling problems. It also identifies 
the analysis and modeling needs for tires. 

The topics covered include: 1) reduction methods for large-scale nonlinear 

analysis, with particular emphasis on treatment of combined loads, displacement- 
dependent and nonconservative loadings; 2) development of simple and efficient 
mixed finite element models for shell analysis, identification of equivalent mixed 

and purely displacement models, and determination of the advantages of using 
mixed models; and 3) effective computational models for large-rotation nonlinear 
problems, based on a total Lagrangian description of the deformation. 


INTRODUCTION 


The problem of tire modeling and analysis has long been an area of major con- 
cern to the tire and aircraft industries. A hierarchy of models varying in the 
degree of sophistication has been proposed. Some of these models are listed in 
Fig. 1 and are sketched in Fig. 2. For a detailed description of the models see 
Ref. 1. The models are grouped into six groups as follows: 

The first group consists of the early tire models which are characterized by 
their simplicity. Among these models are the itttng , beam , and Sling on dasttc (oT 
vi&COeJLcU>tic) foundations. These models were used by Clark and co-workers (Ref. 2). 
Their major drawbacks are: 1) they require extensive experiments to evaluate 
the equivalent properties, and 2) their accuracy and range of validity are not 
known in advance. 

The second group consists of the CO h.d- network models, which are sometimes 
referred to as netting analysis , wherein the inflation pressure is assumed to be 
carried exclusively by the cords (see Ref. 3) . These models have the drawback of 
neglecting both the bending in the tire and the stiffening effect of the rubber. 

The third group of models are the mcmbAane. models , which are based on the use 
of a linear or nonlinear momentless theory of shells (Refs. 4, 5 and 6). Their 
major drawback is that they cannot handle discontinuities in loading, geometry or 
material properties. 

The fourth group is the two-dimensional axis ymmelaic models (Ref. 7), which 
are limited to axisymmetric loadings. 

The fifth group is the thAee- dimensional continuum models. Two approaches 
have been proposed for the analysis of these models. The first approach is based 
on using semi-analytic techniques to reduce the dimensionality of the problem 
(e.g., Fourier expansions in the circumferential direction). The second approach 
is based on using three-dimensional isoparametric solid elements. 

The sixth group of models includes a variety of two -dimensional ihin and thick 
shell models (see, for example. Refs. 8 and 9). Thin shell models neglect trans- 
verse shear deformation, and their use for modeling tires is therefore 
questionable. Anisotropy results in increasing the size of the analysis model, 
and consequently many investigators neglect its effects by using an orthotropic 
model. 

The present paper focuses on the use of two-dimensional thick shell models. 
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DIFFERENT UK NOBELS 


MODEL 

DRAWBACKS 

• SIMPLE MODELS 

STRING 1 ELASTIC 

BEAM ON OR 

RING 1 VISCO- 

ELASTIC 

FOUNDATION 

• REQUIRE EXTENSIVE EXPERIMENTS TO 
EVALUATE EQUIVALENT PROPERTIES 

• ACCURACY AND RANGE OF VALIDITY 
NOT KNOWN IN ADVANCE 

• CORD-NETWORK MODELS 

• NEGLECT BENDING IN TIRE AND 
STIFFENING EFFECTS OF RUBBER 

• MEMBRANE MODELS (BASED ON 
MOMENTLESS THEORY OF SHELLS) 

• CANNOT HANDLE DISCONTINUITIES (OR 
SHARP CHANGES) IN LOADING, GEOMETRY 
(CURVATURE) OR MATERIAL PROPERTIES 

• TWO-DIMENSIONAL AXISYMMETRIC MODELS 

• LIMITED TO AXISYMMETRIC LOADING 

• THREE-DIMENSIONAL MODELS 

• SEMI-ANALYTIC SOLUTIONS 

• THREE-DIMENSIONAL SOLID ELEMENTS 

• COMPUTATIONALLY EXPENSIVE 

• LAMINATED ANISOTROPIC SHELL MODELS 
(THIN AND THICK) 

• THIN SHELL MODELS NOT ADEQUATE 

• EFFECT OF ANISOTROPY CAN BE 
SIGNIFICANT 


Figure 1 



Figure 2 
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TIRE CONFIGURATION AND COMPONENTS 


Typical configurations and components of modern tires are shown in Fig. 3 (see 
Ref. 10). Co mm ercially successful tires are now built as a series of layers of 
flexible high-modulus cords encased in a low-modulus rubber or rubber-like material. 
Hence, a laminated (or layered) model is needed. 
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TIRE LOADS 

The three types of loads applied to the tire and their major characteristics 
are listed in Fig. 4. The three load types are: 

1) In^tcutcon pHQA&USUL, which is axisymmetric but is displacement dependent 

2) Me.ckaru.caZ Zoadi which include centrifugal force, impact loading, contact 
forces, and frictional forces; except for the centrifugal force, which is axisym- 
metric (and displacement dependent) , all the other loads are symmetric 

3) TkcnmaZ ZoacU, which arise due to various manufacturing and operating 
conditions, such as unequal expansion and contraction of rubber and cord, 
hysteretic heating, sliding of the tread on a rough surface, and cord shrinking 
after molding 


LOADS 

CHARACTERISTICS 

• INFLATION PRESSURE 

• AXISYMMETRIC BUT DISPLACEMENT- 
DEPENDENT 

• MECHANICAL LOADS 

• CENTRIFUGAL FORCE 

• IMPACT LOADING 

• CONTACT FORCES 

• FRICTIONAL FORCES 

• AXISYMMETRIC BUT DISPLACEMENT- 
DEPENDENT 

• ASYMMETRIC 

• 

• 

■ 

i 

• THERMAL LOADS 

• UNEQUAL EXPANSION AND 
CONTRACTION OF RUBBER 
AND CHORD 

• HYSTERETIC HEATING 

• CHORD SHRINKING AFTER 
MOLDING 

• ASYMMETRIC 


Figure 4 
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CHARACTERISTICS OF EFFECTIVE SHELL ELEMENTS 
FOR ANALYZING TIRES 

The characteristics of an effective shell finite element model for analyzing 
tires are listed in Fig. 5. The shell element is developed using either a con- 
sistent two-dimensional shell theory or a three-dimensional continuum theory with 
proper interpolation functions in the thickness direction. The elements obtained by 
using the latter approach are referred to as de.ge.ne/iate !>heZZ QZe.m2.nt!>. If a two- 
dimensional shell theory is used, the elements need to be deep and curved and must 
account for each of the following effects: 

1) Laminated construction and anisotropic material behavior 

2) Variation in geometry (e.g., curvature and thickness) as well as of other 

lamination parameters 

3) Transverse shear deformation 

4) Large rotations 

5) Pressure stiffness (for displacement-dependent loadings such as inflation 

pressure) 

6) Thermoviscoelastic material response 


• BASED ON EITHER 

• CONSISTENT TWO-DIMENSIONAL SHELL THEORY, OR 

• THREE-DIMENSIONAL CONTINUUM THEORY WITH PROPER INTERPOLATION FUNCTIONS 
IN THE THICKNESS DIRECTION (DEGENERATE SHELL ELEMENTS) 

• DEEP, CURVED ELEMENTS 

• INCLUDE EFFECTS OF: 

• LAMINATED CONSTRUCTION AND ANISOTROPIC MATERIAL BEHAVIOR 

• VARIATION IN GEOMETRY (E.G,, CURVATURE AND THICKNESS), LAMINATION PARAMETERS 

• TRANSVERSE SHEAR DEFORMATION 

• LARGE ROTATIONS 

• PRESSURE STIFFNESS (FOR DISPLACEMENT-DEPENDENT LOADING; E.G., INFLATION 
PRESSURE) 

• THERMOVISCOELASTIC MATERIAL RESPONSE 


Figure 5 
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OBJECTIVES AND SCOPE 

The objectives of this paper are listed in Fig. 6. They are: 

1) To review some recent developments in finite element technology which are 

applicable to the analysis and modeling of tires 

2) To identify some of the analysis and modeling needs for tires 

The paper is divided into four parts. The first part deals with new develop- 
ments in reduction methods for nonlinear problems. These include computational 
procedures for handling combined, displacement-dependent, and nonconservative 
loads. The second part of the paper deals with mixed finite element models for 
tires in which the fundamental unknowns consist of both force and displacement 
parameters. The equivalence of some of these models with some of the purely 
displacement models is discussed. ' 

The third part of the paper deals with large-rotation nonlinear problems. Two 
formulations are presented; namely, a mixed formulation and a penalty formulation. 
Both formulations are based on the total Lagrangian description of the deformation. 
The fourth and last part of the paper deals with analysis and modeling needs for 
tires. 


OBJECTIVES 

• REVIEW SOME RECENT DEVELOPMENTS IN FINITE ELEMENT TECHNOLOGY WHICH ARE 
APPLICABLE TO ANALYSIS AND MODELING OF TIRES 

• IDENTIFY ANALYSIS AND MODELING NEEDS FOR TIRES 
SCOPE 

• REDUCTION METHODS FOR NONLINEAR PROBLEMS 

• COMBINED LOADING PROBLEMS 

• DISPLACEMENT-DEPENDENT AND NONCONSERVATIVE LOADING PROBLEMS 

• DYNAMIC PROBLEMS 

• MIXED FINITE ELEMENT MODELS 

■ EFFICIENT AND ACCURATE MIXED MODELS 

• EQUIVALENT CLASSES OF MIXED MODELS AND REDUCED/SELECTIVE INTEGRATION 
DISPLACEMENT MODELS 

• MERITS OF MIXED MODELS OVER EQUIVALENT DISPLACEMENT MODELS 

t LARGE-ROTATION NONLINEAR PROBLEMS 

• MIXED FORMULATION 

• PENALTY FORMULATION 


Figure 6 
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REDUCTION METHODS FOR NONLINEAR TIRE PROBLEMS 


The first topic considered in this paper is reduction methods for nonlinear 
analysis. The basic features of reduction methods are outlined in Fig. 7. They 
are techniques for reducing the number of degrees of freedom through the trans- 
formation shown in the figure. The vector {X} represents the original displacement 
degrees of freedom. The vector refers to amplitudes of displacement modes and 
[r] is a transformation matrix whose columns represent a priori chosen global 
displacement modes. 


As is to be expected, the effectiveness of reduction methods depends to a great 
extent on the proper selection of the displacement modes. In a number of studies 
it was shown that an effective choice of the displacement modes includes the various- 
order derivatives of the displacement vector with respect to the load parameter (see 
Refs. 11 and 12). These vectors are generated by uuslng the fantte element model o£ 


2 

3X 3 x 

the tine. The recursion formulas for evaluating the derivatives {-r — } , { — ;r}, ... 

8p 8p 2 

are obtained by successive differentiation of the original finite element equations. 
The left-hand sides of the recursion formulas are the same (see Ref. 12). Therefore 


only one matrix £ actorlzation aj> required far the generation o& all the global 
approximation vectors. Several numerical experiments have demonstrated the effect- 
iveness of this choice (see Refs. 12 and 13). 


DEFINITION: ARE TECHNIQUES FOR REDUCING THE NUMBER OF D.O.F. 
THROUGH THE TRANSFORMATION 

1X, n.l = [r] ^r.l . r < < n 

|X| = ORIGINAL D.O.F. IN THE FINITE ELEMENT 
MODEL 

[H = MATRIX OF GLOBAL DISPLACEMENT MODES 

1*1 = REDUCED D.O.F. -AMPLITUDES OF DISPLACEMENT 
MODES 


JUSTIFICATION: FOR MANY TIRE PROBLEMS THE LARGE NUMBER OF 

D.O.F. |X| IS DICTATED BY THE COMPLEX TOPOLOGY 
OF THE TIRE (DISCONTINUITIES IN GEO/VE TRY, 
LAMINATION, ETC.) RATHER THAN BY EXPECTED 
COMPLEXITY OF BEHAVIOR 


SELECTION OF GLOBAL DISPLACEMENT MODES : 


— 

jdXj jd2x ( 

I 

d 3 X| 

' 


W Lp Z l 

1 

OP 3 [• 


= 

LOAD PARAMETER 



• COLUMNS 0F[r] GENERATED BY USING THE ORIGINAL FINITE ELEMENT 
MODEL OF THE TIRE 

• THEIR GENERATION REQUIRES ONLY ONE 'LARGE MATRIX' FACTORIZATION 

• NUMERICAL EXPERIMENTS HAVE DEMONSTRATED THEI R EFFECTIVENESS 


Figure 7 
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BASIC EQUATIONS USED IN REDUCTION METHODS 
FOR NONLINEAR TIRE PROBLEMS 


The basic equations used in the reduction methods for geometrically nonlinear 
tire problems are given in Fig. 8. It is worth noting that the original displace- 
ment unknowns {X} can be on the order of thousands whereas the reduced unknowns { t|j } 
are typically 20 or less. This is true regardless of the complexity of the structure 
and/or loading. The details of the computational procedure for tracing the load- 
deflection paths in geometrically nonlinear static analysis are given in Refs. 11 
and 12. 



ACTUAL (LARGE) PROBLEM 

REDUCED (SMALL) PROBLEM 

FUNDAMENTAL 

UNKNOWNS 



1X1= INDIVIDUAL - 1X1= [I - ]}*} - AMPLITUDES 

DISPLACEMENTS 
• THOUSANDS OF UNKNOWNS 

Ul- U 1 hr LAUUVltf'll 

MODES 

• TWENTY OR LESS 

GOVERNING 

EQUATIONS 

[K]|Xl + IG (X)|- plPl = 0 
~ 1000 EQUATIONS 

CKW + {G (^)f-plPf = 0 
~ 20 EQUATIONS 

HOW TO TRACE L0AD- 
DEFLECTION PATH 

• REPEATED SOLUTION OF 
LARGE SYSTEMS OF 
SIMULTANEOUS NONLINEAR 
ALGEBRAIC EQUATIONS 

• GENERATION OF [I - ] 

• MARCHING WITH SMALL 
SYSTEM OF EQUATIONS 

• ERROR SENSING AND 
CONTROL (UPDATING [T] 
WHENEVER NEEDED) 


[K] = cr] T CK][n . ig(^)} = [r3 T iG(X)( . ip( = [r] T ipi 


Figure 8 
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APPLICATION OF REDUCTION METHODS TO GEOMETRICALLY NONLINEAR 
ANALYSIS OF A TIRE SUBJECTED TO UNIFORM INTERNAL PRESSURE 

As a simple application of reduction methods to the geometrically nonlinear 
analysis of tires, consider the laminated anisotropic elliptic toroidal shell 
shown in Fig. 9. Due to axial symmetry, only one meridian was modeled using 
four-noded elements with cubic Lagrangian interpolation functions for all the 
displacement and rotation degrees of freedom. The high accuracy of the total 
strain energy obtained by using six basis vectors is demonstrated in Fig. 9. 



G tt = 270 psi h = 0.42 in. 

Fjj = 0.4 e = 0.5 


CORD ORIENTATION + 45/- 45/ + 45/ - 45 ... . 

NL = 10 

Figure 9 
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TREATMENT OF COMBINED LOADS 


The basic equations and the computational procedure used in applying reduction 
methods to the analysis of tires subjected to combined loads are highlighted in 
Figs. 10 and 11. For simplicity, only two independent loads are considered. 

First, the original finite element equations are given. The external loading 
is normalized with respect to two independent load parameters p^ and P 2 ■ The basis 
reduction is done as before, via the transformation shown in Fig. 10. Then the 
Rayleigh-Ritz technique is used to approximate the original set of finite element 
equations by a reduced system of equations in the new unknown parameters {iJj}. The 
number of these equations is considerably less than that of the original equations. 

As previously noted, the crux of reduction methods is the proper selection of 
the transformation matrix [r]. In the case of combined loading, the columns of the 
matrix [f] are selected to be the various-order derivatives of the displacement 
vector {X} with respect to the two independent parameters p^ and 

To trace the different nonlinear paths, corresponding to different combinations 
of the independent load parameters, the basis vectors are evaluated for the unloaded 
structure (p-^ = P 2 = 0) , and the corresponding reduced equations are generated. The 
different nonlinear paths of the tire are obtained by fixing one of the load para- 
meters, varying the other, and repeating the process with different values of the 
first load parameter. This is all done using the -icune. ofa H.Q,dac.^.d ZqucutLoni . 

The total cost of the analysis, to a first approximation, is little more than the 
cost of one linear solution of the original, full system of finite element equations. 
The procedure is described in detail in Ref. 14. 

As a by-product of this technique, a considerable reduction can be made in the 
size of the analysis model used in studying the nonlinear response of tires sub- 
jected to asymmetric loading. This can be accomplished by decomposing the loading 
into symmetric and antisymmetric components and treating each as an independent 
loading. 


REDUCTION METHODS FOR NONLINEAR PROBLEMS 
TREATMENT OF COMDINED LOADS 


GOVERNING FINITE ELEMENT EQUATIONS 


[K){X} + {G (X) } - Pl {p (1) } - p 2 (p (2) } = 0 

[K] = LINEAR GLOBAL STIFFNESS MATRIX 
{X> = VECTOR OF NODAL DISPLACEMENTS 
{G (X) } = VECTOR OF NONLINEAR TERMS 
{P (1> }, {P (2) } = NORMALIZED LOAD VECTORS 
P 1( p 2 = INDEPENDENT LOAD PARAMETERS 

BASIS REDUCTION 




p 


2 


{X} 


n, 1 


m 


n, r ^ } r , ! 


r << n 


REDUCED SYSTEM OF EQUATIONS 

[K]{ljl} + (GW) - Pl {p (1) ) - p 2 {P (2) } = o 


Figure 10 


TREATMENT OF COMBINED LOADS 

SELECTION OF BASIS VECTORS 


(x) = trirti 


in = 


■3P1 


{-52L} 

V> 


3pf 


{• 


} { 


3p l 3p 2 3p 


3 2 X t 

2 / ... 


COMPUTATIONAL PROCEDURE 


• EVALUATE BASIS VECTORS AT p 2 = p 2 = 0 (UNLOADED TIRE) AND GENERATE 
REDUCED EQUATIONS 

• TRACE DIFFERENT EQUILIBRIUM PATH BY FIXING ONE OF THE LOAD 
PARAMETERS AND VARYING THE OTHER (USING T HE SAME SET OF 
REDUCED EQUATIONS) 




p 2 


x 




NOTE ; THIS APPROACH CAN BE USED TO REDUCE THE SIZE 

OF ANALYSIS MODELS FOR THE CASE OF UNSYMMETRIC LOADINGS. 


1 2 


Figure 11 
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TREATMENT OF DISPLACEMENT-DEPENDENT 
AND NONCONSERVATIVE LOADING 

The basic equations used in applying reduction methods to the nonlinear analysis 
of tires subjected to displacement-dependent loading are given in Fig. 12. 

First, the governing finite element equations for the total Lagrangian formu- 
lation are shown. The only new term in these equations is the pressure stiffness 
matrix which represents the follower-load effect, and is unsymmetric for noncon- 
servative loadings. The basis reduction is done and the reduced equations are 
obtained in the manner outlined previously. The following two important facts are 
to be noted: 

1. The basis vectors are evaluated for the unloaded structure. Hence, the 
pressure stiffness matrix does not enter into the left-hand side and only the UneaA 
symmetric global stl^ness matnlx needs to be decomposed. 

2. Since the reduced equations are small in number (on the order of ten or 
less) no symmetrization is needed in the case of nonconservative loadings. 


GOVERNING FINITE ELEMENT EQUATIONS 

FOR A TOTAL LAGRANGIAN FORMULATION 

[[K] - p[K (p) l] {X} + {G (X) } - p{P} = 0 

[K (P) 1 = PRESSURE STIFFNESS MATRIX (UNSYMMETRIC FOR 
NONCONSERVATIVE LOADING) 

BASIS REDUCTION 

(x) = [r] {^} 

REDUCED SYSTEM O F EQUATIONS 

J[K] - P[K (P) 1]U} + (G(<|D) - p(P) = 0 

NOTES : 

• BASIS VECTORS ARE EVALUATED AT P = 0. THEREFORE, ONLY THE SYMMETRIC Ik] 
MATRIX NEEDS TO BE DECOMPOSED. 

• REDUCED EQUATIONS ARE SMALL IN NUMBER (~10). THEREFORE, NO SYMMETRIZATION 
IS NEEDED, 

Figure 12 
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APPLICATION OF REDUCTION METHODS TO THE BIFURCATION BUCKLING 
ANALYSIS OF A RING SUBJECTED TO HYDROSTATIC PRESSURE 


As a simple application of reduction methods to structures subjected to dis- 
placement-dependent loadings, consider the circular ring subjected to hydrostatic 
pressure shown in Fig. 13. 

Doubly-symmetric buckling modes are considered*, hence, only one quadrant of the 
ring was analyzed using higher-order shear-flexible elements with a total of 59 non- 
zero degrees of freedom. The lowest three buckling loads obtained using three, 
four and five vectors are listed in Fig. 13. The lowest buckling load obtained by 
using four vectors agrees, to five significant digits, with that obtained using the 
full system of equations. With five vectors, the error in the third buckling load 
is less than 3%. 


NUMBER OF 
BASIS 

EIGENVALUES p 

P R 3 
K o 

' El 

VECTORS 

p i 

P 2 

P 3 

3 

2. 9998 

15.1888 


4 

2.9997 

14.9964 

36.5126 

5 

2.9997 

14.9957 

35.9668 

FULL SYSTEM 
(59 D.O.F.) 

2.9997 

14.9947 

34.9751 
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REDUCTION METHODS FOR NONLINEAR DYNAMIC PROBLEMS 


The application of reduction methods to transient tire problem is highlighted 
in Fig. 14. First, the governing semi-discrete finite element equations are given 
for the case of no damping. Then the key elements for an effective reduction 
method are listed (see Ref. 13). They include: 

1) The proper selection of basis vectors (the columns of the matrix [r]) 

2) Characterization of nonlinear dynamic response by means of one or few 

scalars 

3) Sensing and controlling the error in the reduced system of equations 


GOVERNING SEM 1- DISCRETE FINITE ELEMENT EQUATIONS 
[Ml IX} = \P\ { - f KJ IX} - {G(X){ 


KEY ELEMENTS FOR EFFECTIVE REDUCTION METHOD 


PROPER SELECTION OF BASIS VECTORS 

{X} n.l =[r, n.r^l' r « " 

CHARACTERIZATION OF NONLINEAR DYNAMIC RESPONSE BY 
MEANS OF ONE OR FEW SCALARS 


• SENSING AND CONTROLLING THE ERROR IN THE REDUCED 
SYSTEM OF EQUATIONS 


Figure 14 
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SELECTION OF BASIS VECTORS FOR THE CASE OF STEP LOADING 


A particular choice of basis vectors which was found to work well for the case 
of step loading is shown in Fig. 15. The vectors consist of a few eigenvectors of 
the linear problem and a few eigenvectors of the steady-state (static) nonlinear 

3G. 

problem. The matrix [ ■ ] is obtained by using the steady-state (static) nonlinear 

OA . 

3 

solution. Reduction methods can be used to reduce the computational effort required 
for generating the steady-state nonlinear solution. 


P(t) 


BASIS VECTORS CONSIST OF: 

• FEW EIGENVECTORS OF LINEAR PROBLEM 

[ K 1 {XI = X [ M 1 {XI TIME t 


• FEW EIGENVECTORS OF STEADY-STATE (STATIC) NONLINEAR PROBLEM 


( K 1 


6G i 


dX. 

J 


{XI 


Xf M ) {X| 


Figure 15 
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APPLICATION TO CLAMPED SHALLOW SPHERICAL CAP SUBJECTED 
TO A CONCENTRATED LOAD AT THE APEX 


As a simple application of reduction methods to nonlinear dynamic problems, 
consider the clamped spherical cap subjected to a concentrated load which has a step 
variation in time (see Fig. 16) . The displacement time history obtained using the 
full system of finite element equations, the reduced system with ten linear vibra- 
tion modes, and the reduced system with the proposed set of modes are shown in 
Fig. 16. The. baAtA ve.ctosu> [ elge.nmodeA ) weJie not updated throughout the. anatyiti . 

As can be seen from Fig. 16, the proposed set of basis vectors predicts qualitatively 
the correct response. The phase shift was almost eliminated by increasing the 
number of basis vectors to 14. 



TIME (MICROSECONDS) 

Figure 16 
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MIXED FINITE ELEMENT MODELS FOR TIRES 


The second topic considered in this paper is mixed finite element models for 
tires. The basic features of the mixed models are outlined in Fig. 17. The finite 
element models include the effects of both laminated anisotropic construction and 
transverse shear deformation, and allow the geometric and material properties to 
vary within individual elements. The fundamental unknowns consist of the eight 
stress resultants and the five generalized displacements. The stress resultants 
are discontinuous at element interfaces, and therefore can be eliminated on the 
element level. 


• TIRE MODELED USING LAMINATED ANISOTROPIC, SHEAR-FLEXIBLE, DEEP SHELL 
ELEMENTS WITH VARIABLE GEOMETRIC AND MATERIAL PROPERTIES 

• FUNDAMENTAL UNKNOWNS ARE-. 

STRESS RESULTANTS N hB , M aB , Q a 

GENERALIZED DISPLACEMENTS u a , w , <f> a 

• STRESS RESULTANTS ARE DISCONTINUOUS AT ELEMENT INTERFACES - ELIMINATED 
ON ELEMENT LEVEL 
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MATHEMATICAL FORMULATION FOR THE MIXED MODEL 

The mathematical formulation for the two-field mixed model is based on the use 
of a moderate-rotation nonlinear shell theory in conjunction with the Hellinger- 
Reissner mixed variational principle. The basic features of this formulation are 
outlined in Fig. 18. Different approximation functions are used for each of the 
stress-resultant fields and the generalized displacement field. The governing finite 
element equations for individual elements can be partitioned as shown in Fig. 18. 

The vector {M} is quadratic in {X} and the vector {G} is bilinear in {H} and {X}. 

For mixed models with discontinuous stress resultants at element interfaces, 
the stress resultants can be eliminated on the element level and the governing finite 
element equations reduce to cubic equations in {X} (see Ref. 15). 

APPROXIMATION FUNCTIO NS 
STRESS RESULTANTS 


= OVllH} . jH|= VECTOR OF STRESS N 

RESULTANT PARAMETERS 

DISPLACEMENTS 


= 01/31X1 . 1X1= VECTOR OF NODAL 

DISPLACEMENTS 

G OVERNIN G FINITE ELEMEN T EQ UATIONS FOR 

■ -cf][s] 1 + n/wx)} vr<n 

,[S]T [OjJlxj + (J<5(H.X)}J W 
WHERE { /Vf ( X ) } AND |G(H, X)| ARE VECTORS OF NONLINEAR TERMS 
(QUADRATIC AND BILINEAR IN jH) AND IX}). 

DISCONTINU OUS STR ES S RESUL TA NTS AT ELEMENT INTERFACES 
1 H } = [F] _1 [S]{X| + [FT 1 \mX)\ 

AND GOVERNING FINITE ELEMENT EQUATIONS REDUCE TO 
[S] T [F] 'I CSJlXf + |G(X)| = |Pf 
WHERE lG(X)f = VECTOR OF NONLINEAR TERMS (CUBIC IN IX}). 


Figure 18 
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REDUCED /SELECTIVE INTEGRATION DISPLACEMENT MODELS 


In recent years a class of displacement models with a performance comparable 
to that of the mixed model has been developed. These are referred to as reduced/ 
selective integration displacement models (see, for example. Refs. 15 and 16). 

The major features of these models are outlined in Fig. 19. The governing finite 
element equations for the individual elements are cubic in {X}. The definitions 
of full, reduced, and selective integration are given in Fig. 19. If an nxn Gauss- 
Legendre formula is used to integrate the linear stiffness matrix [k] exactly for 
parallelogram elements, then in full integration nxn quadrature points are used. 

In reduced integration (n-1) x (n-1) quadrature points are used, and in selective 
integration nxn quadrature points are used for some terms of [k] and (n-1) x (n-1) 
points for other terms. 


GOVERNING FINITE ELEMENT EQUATIONS FOR INDIVIDUAL ELEMENTS 
[K]|X| + {G(XU = {Pi 

IX} = VECTOR OF NODAL DISPLACEMENTS 
|G(X){ = VECTOR OF NONLINEAR TERMS (CUBIC IN |X() 

FULL (NORMAL). REDUCED AND SELECTIVE INTEGRATION 

• IF nxn GAUSS-LEGENDRE FORMULA IS USED TO INTEGRATE [K] EXACTLY 
FOR RECTANGULAR (OR PARALLELOGRAM) ELEMENTS 

• FULL (NORMAL) INTEGRATION USES nxn QUADRATURE POINTS 

• REDUCED INTEGRATION USES (n-1) x (n-1) QUADRATURE POINTS 

• SELECTIVE INTEGRATION USES nxn POINTS FOR SOME TERMS AND 
(n-1) x (n-1) POINTS FOR OTHER TERMS OF [K] 


Figure 19 



EQUIVALENT FINITE ELEMENT MODELS 


The equivalence between finite element models is defined in Fig. 20. Finite 
elements are equivalent if their individual governing equations, when expressed in 
terms of a common set of nodal variables and/or parameters, are identical (see 
Ref. 16). It is important to note that the other parameters not contained in the 
common set are local to the individual elements. For nearly equivalent models, 
the finite element equations are almost identical. 


DEFINITION OF EQUIVALENCE 

• TWO FINITE ELEMENT MODELS ARE EQUIVALENT IF THEIR GOVERNING 
FINITE ELEMENT EQUATIONS, WHEN EXPRESSED IN TERMS OF A COMMON 
SET OF NODAL VARIABLES AND/OR PARAMETERS, ARE IDENTICAL . 

THE OTHER NODAL VARIABLES AND/OR PARAMETERS NOT CONTAINED 
IN THE COMMON SET MUST BE LOCAL TO THE INDIVIDUAL ELEMENTS 
(I.E., DO NOT AFFECT THE ASSEMBLY PROCESS). 

• NEARLY EQUIVALENT MODELS ARE ONES FOR WHICH THE GOVERNING 
FINITE ELEMENT EQUATIONS ARE ALMOST IDENTICAL. 


Figure 20 
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EQUIVALENT MIXED AND DISPLACEMENT MODELS 


The governing finite element equations of both mixed and displacement models 
are shown in Fig. 21 and the mathematical requirements for the equivalence of the 
two models are listed. The table given in Fig. 21 lists examples of equivalent 
quadrilateral (in planform) mixed and displacement models. The following can be 
noted: 

1) Equivalent mixed and displacement models have the same number of displace 
ment nodes, use the same approximation functions for the generalized displacements 
and use the same number of numerical quadrature points. 

2) The symbol (F) refers to full integration and (R) refers to reduced 
integration. 

3) If the geometric and material characteristics within the individual 
elements are constants, the number of quadrature points listed in the table 
generates exact integrals for the mixed models and only approximate integrals 
for the displacement models. 


EQ UIVALENCE 


MIXED MODEL WITH DISCONTINUOUS 
STRESS RESULTANTS 

DISPLACEMENT MODEL 

GOVERNING FINITE ELEMENT EQUATIONS 


CS] T [F3 _1 [S]{X}+{6(X>! = iPi 

[K]{X| + |G(X)| = |P( 

FROM WHICH 


[S] T [F] _1 [S] = 

= [K] 

jS(X)| = 

{G (X)| 


NEAR EQUIVALENCE = IS REPLACED BY - 
EXAMPLES OF EQUIVALENT QUADRILATERAL ELEMENTS 


NUMBER OF 
DISPLACEMENT 
NODES 

QUADRATURE 

POINTS 

NUMBER OF STRESS 
RESULTANT PARAMETERS 


4 

2 x 2 (F) 
1 (R) 

4 

1 


9 

3x 3 (F) 
2 x 2 (R) 

9 

4 


16 

4 x 4 (F) 
3 x 3 (R) 

16 

9 

0 


Figure 21 
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NONLINEAR RESPONSE OF CIRCULAR TOROIDAL SHELL 
SUBJECTED TO UNIFORM EXTERNAL PRESSURE 


To assess the accuracy of the mixed models with discontinuous stress resultants 
at interelement boundaries, the nonlinear response of the circular toroidal shell 
shown in Fig. 22 is analyzed using these models. The solutions obtained using six 
and eight finite elements with nine displacement nodes and four stress nodes are 
compared with the converged solution in Fig. 22. 



E = 1 x 10 7 psl 
p = 100 psl 
a = 15. In. 
b = 10. in. 



x 


Figure 22 
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CLAMPED CYLINDRICAL SHELLS SUBJECTED 
TO UNIFORM PRESSURE LOADING 

To assess the accuracy of the different displacement and mixed models, the 
large-deflection nonlinear response of the clamped cylindrical panel shown in 
Fig. 23 is analyzed using these models. The solutions obtained using a 4x4 grid 
of four-noded quadrilateral elements are shown in Fig. 23. The solutions obtained 
using a 2x2 grid of nine-noded quadrilateral elements are shown in Fig. 24. 

As is to be expected, the full- integration four-noded displacement model is too 
stiff. The full-integration nine-noded displacement model (with the same total 
number of degrees of freedom), though less stiff than the four-noded model, over- 
estimates the stiffness, particularly at higher loads. The mixed model with dis- 
continuous stress resultants is more accurate than the mixed models with continuous 
stress resultants developed in Refs. 17 and 18. 


ACCURACY OF FOUR - NODED QUADRILATERAL ELEMENTS 

4 x 4 GRID 

p o CONVERGED SOLUTION 



Figure 23 
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CLAMPED CYLINDRICAL PANEL 

ACCURACY OF NINE - NODED QUADRILATERAL ELEMENTS 

2 x 2 GRID 

• CONVERGED SOLUTION 



Figure 24 
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USE OF REDUCTION METHODS IN CONJUNCTION 
WITH MIXED MODELS 


The use of reduction methods in conjunction with mixed models is outlined in 
Fig. 25. First, the governing finite element equations for the individual elements 
are given. Then, the vectors of fundamental unknowns (stress resultants and 
displacements) are expressed as linear combinations of a small number of vectors. 
The basis reduction and reduced system of equations are obtained in the manner 
outlined previously. It is important to note that tkz fLQ.du.czd zqucutLonS aAQ 
quadfuxtlc -in the ficduccd unknowns { cj> } . 


GOVERNING FINITE ELEMENT EQUATIONS FOR INDIVIDUAL ELEMENTS 



U} = VECTOR OF UNDETERMINED COEFFICIENTS 
P = LOAD PARAMETER 


REDUCED SYSTEM OF EQUATIONS 


[F] {<(> } + {M (<j>) } = p(P} 


WHERE 


[F] = £ -rr H ] T [F][r H ] + [r H l T [s][r x ] + [ r x ] T [ s ] T [ r R ] 

elements “ L " A 


{M} 


[r H ] T {M(X) } + [r <j ] T {G(H,X) } 


elements 

M>) = VECTOR OF UNDETERMINED COEFFICIENTS OR PARAMETERS 


Figure 25 
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USE OF REDUCTION METHODS IN CONJUNCTION 
WITH DISPLACEMENT MODELS 


The use of reduction methods in conjunction with displacement models, which are 
equivalent to the proposed mixed models, is outlined in Fig. 26. Note that the re- 
sulting reduced system of equations is cubic in the reduced unknowns {<(>}. The 
implication of this is that even l& the mixed model and the dh>plaaement model, axe 
equivalent tkelx xedueed &y&tmi> axe not equivalent . 


GOVERNING FIN ITE ELEMENT EQUATIONS FOR INDIVIDUAL ELEMENTS 
[K] {X} + {G (X)} = pjP} 

BASIS REDUCTION 

M n,l = [ r x] ^r,l * r << n 
n. r 



'a a 1 " 1 

[ r x] = 

lx} dp {X{ . . . . ^ r _j |X{ 


{<&} = VECTOR OF UNDETERMINED COEFFICIENTS 


OR PARAMETERS 
p = LOAD PARAMETER 
REDUCED SYSTEM OF EQUATIONS 
[KJIOI + {G (Q )} = p {P} 

WHERE [K] = [ r X ] T[K] [ r x] 

|G«t>)| = |T x ] t jG OD)} 


= VECTOR OF NONLINEAR TERMS (CUBIC IN (<D)) 


NOTE: EVEN IF MIXED MODEL AND DISPLACEMENT MODEL ARE EQUIVALENT, 
THEIR REDUCED SYSTEMS ARE NOT EQUIVALENT. 


Figure 26 
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ACCURACY OF REDUCTION METHOD - MIXED AND DISPLACEMENT MODELS 

The nonlinear solutions obtained using the reduction method in conjunction 
with the equivalent mixed and reduced- integration displacement models are compared 
in Fig. 27 for the case of a clamped cylindrical shell subjected to uniform pressure 
loading. Seven basis vectors were generated for the unloaded shell. The variations 
of the strain energy with the loading, as predicted by the reduction method 
mixed and displacement models, are shown in Fig. 27. The high accuracy of the 
predictions of the mixed model is clearly seen in this figure. 


FULL SYSTEM (672 D.O.FJ 

▲ MIXED 
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ADVANTAGES OF MIXED MODELS OVER 
EQUIVALENT DISPLACEMENT MODELS 


The advantages of mixed models over the equivalent displacement models are 
listed in Fig. 28. These include: 

1) StmpZtctty ofi fcoAmuZxxtXon. Only quadratic and bilinear terms appear in 
the governing finite element equations. By contrast, the governing equations of 
the displacement model include cubic terms. 

2) If the geometric and material characteristics are constants within each 
element, then mo*t otf the tntdQAaZ* can be. evatuated exactly for the mixed elements 
(even when the element has curved faces and edges) . 

3) The mixed models are bettoA. *vUXed fioA 04 Z UXith Adduction method* in non- 
linear problems, in the sense that: 

a) The basis vectors are simpler to generate. 

b) The mixed models lead to higher accuracy of the solutions obtained by 

the reduced system. This is especially true for stress resultants. 


• SIMPLICITY OF FORMULATION (ONLY QUADRATIC AND BILINEAR TERMS 
APPEAR INGOVERNING FINITE ELEMENT EQUATIONS ) 

• MOST OF THE INTEGRALS CAN BE EVALUATED EXACTLY (EVEN FOR ELEMENTS 
WITH CURVED FACES AND EDGES) 

• BETTER SUITED FOR USE WITH REDUCTION METHODS 

• BASIS VECTORS ARE SIMPLER TO GENERATE 

• BETTER APPROXIMATION PROPERTIES (HIGHER ACCURACY OF 
REDUCED SYSTEM. ESPECIALLY FOR STRESS RESULTANTS) 


Figure 28 
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LARGE ROTATION NONLINEAR PROBLEMS 


The third topic considered in this paper is the large rotation nonlinear 
problems. The basic features of two effective computational models are outlined 
in Fig. 29. In both models a totaZ Lag/uzngta.n description of the deformation is 
used. Consequently, the strain-displacement relations contain tsvigonome&vic 
function* ofi the flotation component* . 

The first computational model is a two-field mixed model with discontinuous 
stress-resultant fields at interelement boundaries. The second model is based on 
the use of the penalty method for handling the trigonometric functions, thereby 
simplifying the analysis. 


FORMULATION 

• TOTAL LAGRANGIAN DESCRIPTION OF DEFORMATION 

• STRAIN - DISPLACEMENT RELATIONS CONTAIN TRIGONOMETRIC FUNCTIONS 
OF ROTATION COMPONENTS 

FINITE ELEMENT MODELING 

• MIXED MODELS WITH DISCONTINUOUS STRESS RESULTANTS 

• PENALTY METHOD FOR HANDLING TRIGONOMETRIC FUNCTIONS 


Figure 29 
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ELASTICA PROBLEM - FORMULATION 


As an application of the proposed computational models, consider the elastica 
problem shown in Fig. 30. In the mixed formulation, the Hellinger-Reissner two- 
field mixed variational principle is used. Transverse shear deformation, though 
small, is included to simplify the formulation. The extensional strain e and the 
transverse shear strain y are trigonometric functions of the rotation <f>. 

The penalty formulation, on the other hand, is based on the Euler-Bernoulli 
type beam theory with both the extensional and transverse shearing strains neglected. 
The axial and transverse displacements u and w are incorporated into the functional 
through the use of constraints and penalty numbers. 


MIXED FORMULATION 


it = /[(Ne+Mk + Qy) 


2 2 2 

1 r N NT cr, . . 

2 EA + eT + GA 


£, y INCLUDE TRIGONOMETRIC FUNCTIONS OF <p 


K _ §± 

ds 



N, M, Q ARE DISCONTINUOUS AT INTERELEMENT BOUNDARIES 


PENALTY FORMULATION 

TT = f Lj El (^ 4 .) 2 + - sin (||) 2 + X 2 - 1 + cos <J>) 2 ] ds 

WHERE X lf X 2 ARE PENALTY NUMBERS 
THE RESULTING STIFFNESS MATRIX IS POSITIVE DEFINITE. 


Figure 30 
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ELASTICA PROBLEM - NUMERICAL RESULTS 


The displacements, rotations and total strain energy obtained by using fiouA. 
two-nodzd. zZtirnnti and Xwo tkA<lZ-nodQ.d eZmwtA are depicted in Figs. 31 and 32. 
Also, the deformed configurations of the beam for various values of the transverse 
load P are shown in Fig. 31. Both exact-integration displacement models (DE models) 
and mixed models with discontinuous forces (MD models) are used. As to be expected, 
the displacement models are too stiff. This is particularly true for the two-noded 
elements. By contrast, the predictions of the mixed models are highly accurate. 


LARGE-ROTATION ELASTICA PROBLEM 



— EXACT 

SOLUTION 

o 

A DE2 ) 

DISPLACEMENT MODELS 

0 

2 DE3 ) 


+ 

A MD2-1 1 

MIXED MODELS 

A 

2 MD3-2 ) 




Figure 31 
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ANALYSIS AND MODELING NEEDS FOR TIRES 


The fourth part of the paper deals with future analysis and modeling needs for 
tires. The overall goal is to develop a general tire analysis capability which 
includes (see Fig. 33): 

1) Accurate representation of the tire configuration and construction 

2) Reliable material characterization including thermo-viscoelastic response 

3) Capability for predicting the stresses and deformations due to footprint 
loading; this also includes the prediction of the contact area 

Since there is a certain degree of uncertainty in the accuracy of the various 
elements of the tire model, considerable work should be directed towards assessing 
the sensitivity of the tire response to various modeling details such as material 
characteristics, surface inaccuracies, and variations in the tire design variables. 
The result of such sensitivity study would allow the identification of the minimum 
degree of sophistication of the model required to achieve a prescibed level of 
accuracy. 

There is also a need for identifying failure mechanisms and developing a 
verifiable failure analysis capability for tires. Use can be made of the consider- 
able experience gained in damage tolerance design concepts for fibrous composite 
structures. 

• GENERAL TIRE ANALYSIS CAPABILITY 

• ACCURATE REPRESENTATION OF TIRE CONFIGURATION AND CONSTRUCTION 

• RELIABLE MATERIAL CHARACTERIZATION INCLUDING THERMO-VISCOELASTIC RESPONSE 

• PREDICTION OF STRESSES AND DEFORMATIONS DUE TO FOOTPRINT LOADING 

(CONTACT AREA, STRESSES AND SLIP) 

• SENSITIVITY ANALYSIS - SENSITIVITY OF RESPONSE TO: 

• MATERIAL CHARACTERISTICS 

• SURFACE INACCURACIES 

• VARIATIONS IN DESIGN VARIABLES (REQUIRED FOR EVALUATION OF STRUCTURAL 

CONCEPTS AND FOR OPTIMIZATION) 

• MODELING DETAILS (IN ORDER TO DEVELOP SIMPLE TIRE MODELS) 

• FAILURE MECHANISMS AND FAILURE ANALYSIS OF TIRES 

• DAMAGE TOLERANCE 


Figure 33 
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SUMMARY 


In summary, four topics are covered in this paper; namely, recent advances in 
reduction methods for nonlinear problems, mixed models for tires, computational 
models for large-rotation nonlinear problems, and analysis and modeling needs for 
tires. (See Fig. 34.) 

Reduction methods have proven to be very effective for the nonlinear static 
analysis of structures subjected to either combined loads or displacement- 
dependent loads. However, more work is needed to realize their full potential for 
nonlinear dynamic and time-dependent problems. 

Mixed shell models with discontinuous stress resultants at element interfaces 
have high potential for nonlinear analysis of tires. These models can be easily 
incorporated into existing general-purpose finite element programs based on the 
displacement formulation. 

Two computational models are presented for the large-rotation nonlinear 
problems. Both models use a total Lagrangian description of the deformation. 

The first model uses a mixed formulation, and the second model uses a penalty 
formulation. Both models appear to have high potential. 

As far as analysis and modeling needs are concerned, three areas have been 
identified. As is to be expected, the modeling and analysis of tires will be strongly 
impacted by new advances in materials technology, computer hardware, software, 
integrated analysis, and CAD/CAM systems. 


• REDUCTION METHODS 

• VERIFIED FOR STATIC NONLINEAR PROBLEMS INCLUDING CASES OF COMBINED 
LOADS AND DISPLACEMENT-DEPENDENT LOADS 

• FURTHER DEVELOPMENT NEEDED FOR NONLINEAR DYNAMIC PROBLEMS 

• MIXED MODELS WITH DISCONTINUOUS STRESS RESULTANTS 

• HAVE HIGH POTENTIAL FOR ANALYZING TIRES 

• CAN BE EASILY INCORPORATED INTO EXISTING GENERAL-PURPOSE FINITE ELEMENT 
PROGRAMS 

• LARGE ROTATION NONLINEAR PROBLEMS 

• BOTH MIXED AND PENALTY FORMULATIONS PROVIDE EFFECTIVE ANALYSIS 
TECHNIQUES 

• ANALYSIS AND MODELING NEEDS 

• GENERAL ANALYSIS CAPABILITY FOR TIRES 

• SENSITIVITY ANALYSIS 

• FAILURE MECHANISMS AND FAILURE ANALYSIS 

Figure 34 
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